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Abstract 

We extend Onsager’s minimum dissipation principle to stationary states that are only subject to 
local equilibrium constraints, even when the transport coefficients depend on the thermodynamic 
forces. Crucial to this generalization is a decomposition of the thermodynamic forces into those 
that are held fixed by the boundary conditions, and the subspace which is orthogonal with respect 
to the metric defined by the transport coefficients. We are then able to apply Onsager and 
Machlup’s proof to the second set of forces. As an example we consider two-dimensional nonlinear 
diffusion coupled to two reservoirs at different temperatures. Our extension differs from that of 
Bertini, et al. in that we assume microscopic irreversibility and we allow a nonlinear dependence 
of the fluxes on the forces. 
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I. INTRODUCTION. RELAXATION TO EQUILIBRIUM 


Equilibrium statistical thermodynamics is able to estimate the macroscopic quantities 
and their fluctuations without solving the equation of motion. However, an analysis of dy¬ 
namical properties of the systems is essential if we are interested in determining the macro¬ 
scopic behaviour of the systems relaxing to non-equilibrium steady states or in obtaining 
the probability distributions of fluctuations. In such cases, the Gibbs distribution has to be 
replaced by a suitable distribution of fluctuations valid for the microscopic dynamics. This 
theory should be based on the extension of the Boltzmann-Einstein equilibrium fluctuation 
theory combined with dynamics. One example of such example is the macroscopic fluctu¬ 
ation theory of Ref. [I]. This theory has been used for studying some microscopic models 
and it leads to various interesting predictions 0-0 In this paper we develop a theory for 
non-equilibrium steady states mainly based on the Boltzmann-Einstein theory. Crucial in 
our approach is the decomposition of the thermodynamic forces that are held fixed by the 
boundary conditions (fixed thermodynamic forces) into those, which have no external inter¬ 
actions (free thermodynamic forces). We shall show that, without explicit knowledge of the 
entire invariant distribution function for the microscopic dynamics, the Onsager-Machlup 
functional, restricted only to the free thermodynamic forces, approximates the probability 
of a particular relaxation to a stationary state. We also provide an estimation of the error. 
The Boltzmann-Einstein theory received a rigorous mathematical formulation in classical 
equilibrium statistical mechanics via the so called large deviation theory (LDT) [B]. The 
LDT has also been applied to hydrodynamic evolutions of stochastic interacting particle 
systems [7j and extended to nonlinear hydrodynamic regime |H]. A general theory of large 
deviations for irreversible processes, i.e. when the detailed balance does not hold, has been 
successively formulated in 2002 by Bertini et al. i-IU. Several examples of LDT are 
provided by thermodynamic systems driven towards non-equilibrium steady-states by the 
boundary conditions. This is the case, for example, of a fluid in contact within two thermal 
reservoirs where a flow of matter, or heat, through the system is established. Bertini et al. 
showed that the spontaneous fluctuations of the process is described by the time reversed 
dynamics [9]. In Ref. |9] it is shown that the violation of the Onsager-Machlup symmetry 
observed, for example, in stochastically perturbed reversible electronic devices [12], is also 
connected to the time reversed dynamics. 
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When the thermodynamic forces, indicated with X'^, are nonzero, entropy s is produced 
according to the balance equation 


ds 

m 


+ V-J 


djs 

dt 


( 1 ) 


where / and djs are the reversible entropy flow and the density of entropy production, respec¬ 
tively. Let us consider now a system characterized by n degrees of variables Ai, A 2 , ■■■ An. 
The (local) equilibrium values are A^, Denoting hy = A^^ — A^^ {i = 1 ■ ■ ■ n) the 

n deviations of the thermodynamic quantities from their equilibrium value (fluctuations), 
Prigogine proposed that the probability distribution of hnding a state in which the values 
lie between and -|- da^j, is given, up to a normalization constant, by [I3] 

^ ocexp(-M) (2) 

with ks denoting Boltzmann’s constant. The negative sign in Eq. ([^ is due to the fact that, 
during the processes, —A/S' < 0. Indeed, if —A/S' were positive, the transformation a* —>■ a' 
would be a spontaneous irreversible change and thus be incompatible with the assumption 
that the initial state is a (local) equilibrium state [13]. The Prigogine theory generalizes 
the Boltzmann-Einstein theory [T3| , [T3] , which applies only to equilibrium thermodynamic 
fluctuations or to adiabatic transformations na, ng, na. The thermodynamic forces, 
and the conjugate fluxes, J^, are related to thermodynamic fluctuations by 

= = (3) 

where the dot stands for the substantial time derivative. Here we have followed Onsager 
and Machlup's notation, treating the forces as if they were scalars as in the case of chemical 
reactions. However the same arguments apply to vector and tensor forces obtained by 
considering as thermodynamic forces the single components of these quantities. 

An equilibrium state is stable if a perturbation in the densities leads to a restoring force 
X^. The restoring force creates a current which cancels the perturbation. Near 

equilibrium Onsager [IB] has extended a theorem of Lord Rayleigh [H] that this process 
is completely characterized by the minimization of a dissipation functional with respect to 
a variation of the currents. In this note we will extend Onsager and Machlup’s stochastic 
derivation of this result 1201 to more general processes in which the system is connected 
to multiple reservoirs at which certain intensive variables are held hxed. These reservoirs 
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prevent the system from arriving at equilibrium, as there will be a flux of the corresponding 
extensive quantities through the system from one reservoir to another. However the system 
will nonetheless, under certain conditions, arrive at a stationary state which is locally in 
equilibrium, and we will extend Onsager and Machlup’s argument to demonstrate that this 
relaxation process is characterized by the minimization of only those parts of the dissipation 
functional corresponding to an orthogonal subset of the forces to the residual current. To 
make this story quantitative, one needs to know the forces as a function of the extensive 
variables and the currents as a function of the forces X'^. While the system is no 
longer in equilibrium, we will assume that locally it is still in equilibrium, and so one may 
locally dehne the total entropy density s. In general, the thermodynamic fluxes are functions 
of the forces. At the thermodynamical equilibrium the thermodynamic forces are zero and 
the fluxes vanish too. The transport equations {i.e., the flux-force relation) may be brought 
into the form 

J, = g,AX)X^ (4) 

where X^ denotes the thermodynamic forces per unit volume and the elements of the matrix 
gfj,u are identihed with the transport coefficients. In Eqs (|^, as well as in the sequel, 
Einstein’s summation convention on the repeated indexes is implicitly understood. Given 
the function g^,^ and a conhguration one may calculate and so determine the entire 
future evolution of the system. When the forces are small, corresponding to small deviations 
from equilibrium, one may approximate gf^^ to be independent of X^, corresponding to fluxes 
Jfj_ which are linear in X^. We refer to this region as the linear Onsager region. In this region, 
Onsager noted that the phenomenological linear relation (|^, which gives as a function 
of maximizes the functional 

M = I s d:s. — - j g^'"a^Ciyd'x. (5) 

Jn 2 

varied with respect to everywhere except for the boundaries (with dx denoting a spatial 
volume element, and the integration is over the entire space G occupied by the system). In a 
dynamical context, one may ask what is the most probable trajectory followed by the system 
in the spontaneous emergence of a fluctuation or during the relaxation towards equilibrium. 
Under the assumption of time reversibility and by using a stochastic argument, Onsager 
and Machlup, demonstrate that the most probable trajectory is obtained by minimizing the 
quantity M in the situation of a linear macroscopic equation (i.e., very close to equilibrium) 
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[20] , [21]. One may now extend the question by asking, for example, what is the most prob¬ 
able trajectory followed by the system in the relaxation to the boundary driven stationary 
non-equilibrium states. In their paper, Bertini et Al. formulated a dynamical fluctuation 
theory for non-equilibrium steady-states, which is based on the time reversed dynamics [9] . 
In the present work we shall modify the Onsager-Machlup theory for stationary non- 
equilibrium states to include cases with a nonlinear macroscopic equation. We shall gen¬ 
eralize the minimum dissipation principle to relaxation processes to steady states that are 
only locally in equilibrium. The hypothesis of local equilibrium guarantees the existence of 
the local entropy functional, as well as the reversibility of the microscopic physics. However 
we do not assume that the transport coefficient gfj^As are independent of the forces and 
so we are not restricted to the linear Onsager region. Our systems will be prevented from 
reaching equilibrium by being coupled to several reservoirs held at distinct values of the 
intensive variables, so that the corresponding forces will be non vanishing in our system. 
These forces will yield fluxes via Q, which will transport some of the quantities between 
the reservoirs. The steady-state satishes the relation satishes the relation 

J^X>^ = = 0 ( 6 ) 

Hence, in a stationary state, some of the currents J^ = a^ = will vanish, those orthogo¬ 
nal to the non-vanishing forces using the metric evaluated at the stationary value of the 
forces. We shall show that if only this second class of quantities = (3^ are perturbed from 
their stationary state, then their relaxation will extremize M varied over only this second 
class of currents. Note that the nonlinear dynamics is nontrivial, because the subspace of the 
currents which is varied is the orthogonal compliment of the forces between the reservoirs, 
and this compliment depends on the full force-dependent metric g^^. The extension of the 
minimum dissipation principle to the nonlinear regime is part of the generalized thermody¬ 
namical program reviewed in Ref. [22]. In particular, a geometric analogue, known as the 
minimum rate of dissipation principle, was introduced in Ref. [23] • A different generaliza¬ 
tion of the minimum dissipation principle to non-equilibrium steady states has appeared in 
Ref. [21]. They continue to work in the linear region, but they relax the local equilibrium 
condition so that the microscopic dynamics is not necessarily reversible. This leads to a 
correction to the flux terms in the quantity M, subtracting out the irreversible part. In¬ 
stead of basing their proof on a stochastic process, which would be difficult away from local 
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equilibrium, they use the Hamilton-Jacobi equations of Ref. i. ra which describe the most 
probable trajectory by which a fluctuation is created. 

In Sec. |IV| we show that even in the nonlinear case the minimum dissipation principle 
determines the evolution of the quantities that relax to equilibrium, and then we adapt the 
Onsager-Machlup derivation of the minimum dissipation principle to the nonlinear case. To 
avoid misunderstandings, the Minimum Dissipation Principle herein expressed should be 
understood as an approximate variational principle in nonlinear transport for systems out 
of equilibrium and not as a variational principle, which is rigorously satished for systems out 
of the linear (Onsager) region. A discussion about an estimation of the error can be found 
in the sub-Sec. IV B Then in Sec. |V]an example is provided which is a nonlinear version of 
that of Ref. [9]. We shall consider a 2-dimensional, boundary driven, nonlinear zero-range 
diffusion in which there are only reservoirs on two opposite faces. We shall see that the 
minimum dissipation principle is obeyed for diffusion along a direction which depends on 
the full nonlinear metric. 

• Equilibrium and Boltzmann’s Principle 

Consider a closed system which is macroscopically characterized by a vector of extrinsic 
quantities Motivated by Boltzmann’s kinetic theory of gases [H], Planck has dehned 
the corresponding entropy S{a) to be, up to a constant shift, proportional to the logarithm 
of the number of microstates in the microcanonical ensemble which yield the macrostate 
described by [26]. In a classical system the number of microstates is inhnite, but he 
suggests that they be coarse grained into discrete quantities which can then be counted. 
This has the advantage that it allows him to accurately describe black body radiation, but 
also that the discretization procedure allows him to dehne a notion of probability (and in 
particular a measure) in an apparently deterministic system. More precisely, he recasts 
Boltzmann’s principle as the ergodicity assumption that, for some choice of discretization, 
the probability of each allowed microstate will be equal. 

In a very large system, to which one may apply the thermodynamic approximation, there 
will be a small region in the space of values of a, peaked about some value for which the 
entropy is much larger than even the integral over the rest of a-space. As is a maximum 
of the entropy function, it satishes 

dAS{a 


da. 


( 0 )= 0 . 


(7) 
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Boltzmann’s principle then states that, given no other information about the system, in the 
thermodynamic limit of a large number of microstates the system will certainly be found in 
the macrostate which is referred to as the equilibrium state. 

• Relaxation in the Onsager Region 

If on the other hand one begins with an initial condition at time t = 0, then 

the probability that the state is in a given state a/ at time t is given by a conditional 
probability, in contrast with the absolute probability in equilibrium. However, due to the 
ergodicity assumption, if one waits an inhnite amount of time then the information about 
the original state is erased and so again the measured state will be the equilibrium 
state with probability 1. If the conditional probabilities can be calculated, then one can 
determine not only the hnal state, but also the full time-dependent trajectory followed by 
the macrostate from to 

How does such a calculation proceed? The tendency of the system to seek equilibrium is 
measured by the thermodynamic forces 




dAiS{a) 

daf. 


0 . 


( 8 ) 


where we have introduced the entropy production A/S', and the fact that at the 

system is not at the local equilibrium state. The fluxes (of matter, heat, electricity) are 
measured by the time derivatives of the a’s: J^(q!) = As in general is a vector, so 
is X^(a). depends on such that when vanishes, so does = J^. The essential 
physical assumption about the irreversible processes is that the fluxes depend on the forces 
through the transport coefficients g^^{X{a),a): 


J^{a) 


da^ 

dt 


g^,{X{a),a)X-{a). 


(9) 


where, here and in the sequel, the summation convention on the repeated indices is un¬ 
derstood. The conditional probabilities can be used to calculate the transport coefficients 
which then reduce the relaxation of the system to the solution the a system of coupled, 
nonlinear, hrst order differential equations (|^. The entropy production is a function of the 
a’s: A/S' = A/S'(ai, • • ■ , a^). Hence, 

dXjS _ dAiS da^ _ 

dt da^ dt ^ 
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In general it is difficult or impossible to determine these conditional probabilities, and so 
this approach is of limited use. However near thermodynamic equilibrium, the entropy 
production is near its extreme value, and it can be expanded about 


/\iS{a) ~ A/5o - = -2s^''a^ (11) 


where is a constant matrix. In this case, using the fact that = 0, one may 

approximate the current to be 

^ a)X''{a) ~ (12) 

This is a system of linear differential equations depending upon two constant matrices. 

• Onsager Machlup Principle 

In Ref. [IH] Onsager has generalized a theorem of Lord Rayleigh [12], that relaxation to 
equilibrium can be derived from a variational principle. They found that Eq. S can be 
derived as an extremum of the quantity with respect to the currents 

M = X^J^ - (13) 


where g^'^ is the inverse matrix of the transport coefficients. We get 


0 = ^ = (14) 

More than 20 years later an interpretation of this formula, near thermodynamic equilibrium, 
was demonstrated. In Ref. m Onsager and Machlup provided a stochastic demonstration 
that close to equilibrium, the probability / of relaxing from ai to Oq along a path a(t) is 
proportional to the exponential of M or more precisely 


/ oc exp 


2 k 


B 


dt {X^J^ - 


. ( 16 ) 


Thus, just as Boltzmann’s principle provides a formula for the absolute probability of re¬ 
alizing a certain microstate at a hxed time, Onsager and Machlup found a formula for the 
probability of following a succession of events during relaxation to equilibrium. 


How did Onsager and Machlup demonstrate that (15) indeed provides the probability of 


any given sequence of conhgurations? Their demonstration rested upon three assumptions. 
First, the sequence of events is Markovian. This means that given some set of data at time 






which in Ref. [20] is the vector a but the entire approach can be generalized to include their 
time derivatives in systems with inertia |2I] , the conditional probability of a configuration at 
any future time ti is independent of any knowledge of the state at times before to- Thus while 
the theory remains nondeterministic because it the configuration at to does not determine 
that at fi, nonetheless it does not contain any hidden variables at to which might affect ti. 
In this sense the a provide a complete set of states. The Markovian property implies that 
to determine the probability of a trajectory, one only needs the conditional probability of 
one state given another. 

The second assumption is that the entropy production function is quadratic, as in Eq. (0. 
This is only a technical assumption to allow a calculation of a path integral in closed form. 
Given a form of the entropy which is not quadratic, perturbation theory could be applied 
to the path integral to yield a result as an asymptotic series, as was already known at the 
time of Onsager and Machlup’s paper. 

Finally they assumed that each extrinsic variable is a sum of local, uncorrelated variables. 
As a result, steps in the evolution of these variables, while random, will obey a Gaussian 
distribution in the thermodynamic limit. This Gaussian distribution is evident in the fact 


that (15) is expressed as an exponent of squares. 


II. RELAXATION TO A STEADY STATE 

For more than half a century there have been attempts to generalize Onsager and 
Machlup’s variational principle to the process of relaxation not to equilibrium, but to a sta¬ 
tionary state. A stationary state is a state satisfying the condition = 0. Hence, at the 

steady state, not all the thermodynamic forces, vanish. Glansdorff and Prigogine have 
argued that relaxation to a stationary state is not described by an extremization problem in 
Ref. [T7]. Indeed in Ref. [27] Silhavy states that in a stationary state the thermodynamic 
quantities themselves are ambiguous. 

The basic problem is the microscopic state is no longer simply an unknown element of 
the microcanonical ensemble. Information about the microstate is constantly destroyed by 
interaction with the outside of the system. Thus instead of being characterized only by the 
internal degrees of freedom, the evolution of a microscopic state also depends on external 
degrees of freedom. Furthermore some of the information about the internal degrees of 
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freedom cease to be measurable on the inside, as they are transfered to the outside, where 
no measurements occur. Thus the starting point of this analysis, a dehnition of the entropy 
as a sum of the number of microstates, is already ill-dehned. 

In this note we will claim that nonetheless a subset of the system’s information is approxi¬ 
mately determined by an extremization principle. We will provide a criterion which describes 
when this approximation is reliable. 


III. EXAMPLE: A PARTIALLY IONIZED PLASMA 

Consider a closed system containing weakly ionized hydrogen gas in a sealed container 
with perfectly reflective walls. The system contains neutral hydrogen atoms, W = W 
free electrons and protons as well as a number of free photons, which we will call iV^, making 
the crude approximation that the energy of each photon is equal to the binding energy of 
an electron in hydrogen. Similarly we will assert that all of the bound electrons are in a Is 
orbital. Collisions and photon absorption cause the hydrogen to ionize, but also electrons 
and protons and recombine into a hydrogen atom and a photon. This photon is just at 
the right energy to ionize another hydrogen, although usually it needs to bounce off of the 
perfectly reflecting walls a few times first. 

The microstates correspond to the discretized positions and velocities of the various particles. 
Summing over these positions one can calculate the independent macroscopic variables Nh, 
Ne and N^. These are components of the vector a. Given some information, like the size 
of the system and the total energy, one can in principle calculate the number of microstates 
corresponding to each value of a, and so determine the functional form of the entropy. 
Maximizing this entropy one can hnd the state which yields, for example, the ionization 
fraction of the hydrogen in equilibrium. It will be approximately given by the Saha equation. 
A short burst with a laser can ionize some atomic hydrogen, taking the system out of the 
old equilibrium. A slight expansion of the cavity can remove this additional energy from 
the system. Now W will be higher than its equilibrium value and W will be lower. One 
can understand the relaxation of this excited state to the new equilibrium using the analysis 
of Onsager and Machlup described above. The excess of free protons and free electrons 
corresponds to a nonvanishing gradient of the entropy function, which is a thermodynamic 
force. This force causes a current, which in this case is just an increase in Nh at the expense 
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of Ni. The total energy is conserved, this change occurs because there are simply more 
states available with a at its equilibrium level Oq than at the excited level. 

How can a nonequilibrium steady state be constructed in this example? Imagine that the 
laser is permanently turned on, but that the plasma is allowed to transfer some kinetic 
energy to the walls of the container. For simplicity we will assume that these collisions are 
sufficiently elastic that they do not affect the rates of ionization and recombination. As 
a result the time derivative of the number of photons will not be determined entirely 
by the thermodynamic forces, there will also be a contribution from the external current. 
However, after waiting for a sufficiently long time, the plasma will heat to a temperature at 
which a = 0. 

What can Onsager’s formalism tell us about this situation? Certainly it cannot tell us 
photon current Nh, which has an external contribution. However the the numbers Nh and 
Ng of hydrogen atoms and free electrons are determined entirely by the internal physics. 
This consists of the same equations as in the case of a closed system, as the recombination 
and ionization processes are unaffected. All that has changed is the number of photons N-y. 
But near the steady state, lies close to a known value. 

In particular the number of classical states available at a given moment S{a) is still well- 
dehned. These are not Hamiltonian eigenstates of a quantum theory which require an 
inhnite amount of time to be manifested, they are coarse grained positions and velocities 
which exist instantaneously. The ergodicity assumption means that the interactions, in this 
case recombination and ionization, will act so as to increase this number of states S. After 
all the rate of ionization and recombination are independent of whether the photon arrived 
from an external source, or was always there bouncing off of the walls. Thus not only the 
energy, but also the thermodynamic forces X and the transport coefficients are given by the 
same formulas as in the case of relaxation to equilibrium. 

The only step in the logical sequence which differs is then the current itself. The time 
derivative of has an external contribution. But the other currents are unaffected, and so 
continue to be given by Eq. (g. 
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IV. THE GENERALIZATION 


A. A Functional for the Internal Quantities 


Now what happens if we perturb the system from a stationary state and try to follow 
Onsager and Machlup’s stochastic derivation of the probabilities of various relaxation pro¬ 
cesses? 

The hrst problem that arises is that in a stationary state, ASj is not extremized. By 
dehnition the time derivatives vanish at a stationary state. However this no longer 
implies that the forces vanish. We can project the vector into two vectors and 7 ^. 
Here 7 ^ represents the quantities which are subject to external sources, and as in a closed 
system, are absolutely conserved. In the plasma example, contains Nh and W and 7 ^ 
contains N^. It is not difficult to show that for any given conhguration we can choose 
a basis of the 7 ^ orthogonal to all of the with respect to the transport coefficients 
Indeed, the evolution equation for fluctuations may be brought into the form 


cLg.^ 

dt 




dAiS 

day 






dAiS 

day 




(16) 


with The greek indexes fi, u, k, ■■■ run from 1, ■ ■ ■ ,n, and the 

second term on the right-hand side of Eq. (16), takes into account the contribution 

to the dynamics due to thermodynamic forces held hxed by the boundary conditions. Note 
that the first m components of the column vector correspond to the m (with m < 

n) thermodynamic forces held hxed by the boundary conditions, whereas the remaining 
n — m components of this vector are set to zero. The vector a^ may be projected into two 
huctuating vectors (3^ and 7 ^ by introducing two matrices and which are orthogonal 
with respect to the metric dehned by the transport coefficients. In particular, we dehne 


X = A’'a. 


IfJ. — 


(17) 


with matrices and satisfying the orthogonal conditions 


.47„B:! = 0 and 


(18) 


The matrices A^ and are related to their inverse matrices through the relations 


g^^A'^^g^^Al = A>^^A^y = and g'^^B^^g^^Bl = B>^^B^y = 




(19) 
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with 6 ^ denoting Kronecker’s delta and the metric tensor is used for raising and/or 
lowering the indexes. Note that from Eq. (19) we get the useful relations 


and B g^p 


( 20 ) 


By taking into account Eqs (17) and (20), the symmetry relation g^p = g^^j,, and the identity 

( 21 ) 


dAjS ^^^dAjS ^ ^^dAjS 


da,, 


d[5p 


dlv 


it is easily checked that in the basis (3^ and 7 ^, the evolution equations (16) factorize 

dAjS 


df3^ 

dt d(3p, 


d^ _ 

dt 


_l_ nu j(ext.) 


( 22 ) 


As the evolution of 7 ^ is affected by external forces, it cannot be determined even statistically 
by a thermodynamic argument of the system itself HU- The most naive approach would be 
to simply apply Onsager’s formalism to /3^, ignoring the 7 ^. After all the evolution of in 
a relaxation to a steady state is described by equations which look like those describing the 
evolution of a^ in a relaxation to equilibrium. 

Recall that Onsager and Machlup’s derivation rested upon three pillars. The third, that the 
extrinsic variables are sums of uncorrelated variables, is unaffected. The dehnitions of these 
variables do not depend upon the external forces. 

However there is one critical difference between the equations describing and In the 
later case, the entropy production A/S'(/ 3 ^, 7 ^) is a function not just of but also of 7 ^. 

This is a problem for the hrst two pillars. The are no longer a complete set, and so the 
Erst pillar, the Markovian property featured in Onsager and Machlup’s derivation is lost. 
The 7 ^ are hidden variables and they create a systematic bias not only in the statistical 
fluctuations, but also in the mean evolution of the (3^. What about the second pillar? The 
entropy still extremizes and in a stationary solution 7 ^ is fixed to some value 7 °. Setting 
7 ^ = 7 °, the action may still be expanded in (3^ and it will be quadratic. However 7 will 
change during the relaxation process, and likely during the excitation away from the steady 
state. Thus there will be corrections to the entropy which depend upon 7 ^ — 7 °. 

Both of these problems have the same cause. The evolution of 7 ^ during the relaxation is 
unknown. It depends not only on the internal physics, but also on the external source. As 
the entropy couples {3^ and 7 ^, this causes a hnite error in an estimate of the trajectory of 
the /?^’s alone using Onsager and Machlup’s functional. 
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B. Estimating the Error 


Just how big is this error? When must it be considered? Of course it depends on just 
how much 7 deviates from its stationary value 70, which depends on a physical choice. We 
are interested in a system which is excited from a steady state and then decays. In the case 
of our partially ionized plasma subjected to a constant laser, for example the laser may be 
turned off for a moment, taking the system out of its steady state. There will momentarily 
be some additional recombination and the temperature will drop. But eventually, as less 
kinetic energy will be transferred to the walls, the system will relax to its steady state. 

In this case the dependence upon the unknown physics of the external photon injection is 
essential. The perturbation from the steady state itself changes the photon number N^, 
which is 7 in this example. It is essential to know just how 7 returns to 70. This restoration 
depends strongly on the physics of the photon injection into the system. It cannot, to any 
approximation, be described by the internal physics of the system. Therefore, an application 
of the Onsager Machlup to the evolution in this case appears to be at best difficult, as the 
information from the functional needs to be supplemented by information about the laser. 
On the other hand, one can also consider a perturbation from a steady state in which one 
only changes /?, leaving 7 = 70. If 7 continues to be equal to 70 during the entire relaxation 
process, then the Onsager Machlup functional applied to (3 will function as well as in a 
relaxation to equilibrium. Of course 7 will not be equal to 70 throughout the relaxation, 
even if it is equal at the beginning and the end. It is this transitory deviation of 7 from 70 
which contributes an unavoidable error to the Onsager Machlup functional for the probability 
of a given relaxation to a steady state. 

We can estimate the error by expanding the entropy to leading order about the steady state 
values of (3 and 7, which for simplicity we will set to /do = 7o = 0. The entropy production 

is 

A/S' = I3al3 + 76/? + 7C7. (23) 

Now we will perturb /3 to /3i. This will exert a force 

A' = ?N^ = ()ft+2c7. (24) 

(77 

The force vanishes when 

7 = -^c“^6/?i. (25) 
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However 7 is also effected by an external force, which compensates for this effect and so 


generically can cause an appreciable shift in 7 from the value in Eq. (25). However, while 


this shift can easily be smaller than or comparable to the value in (25) itself, without hue 


tuning it will not be much greater. Therefore Eq. (25) provides a rough upper limit on the 


7 which can be expected during the relaxation. This will be sufficient for our estimation. 
Now that during the relaxation 7 will shift, so will the action for f3 itself. This backreaction 
of 7 upon /3 is the obstruction to a variational principle describing the relaxation to a steady 


state. To estimate it, we can simply insert Eq. (25) into Eq. (23) 


ArS = ArS^ + 6 ArS 


AjS^ = Pal3 ; 5AiS = (3hc-^bl3 


(26) 


Here A^S^ is the contribution which, using the Stochastic derivation of Ref. [20], would 
lead a probability for each relaxation. The term 6 AjS is the correction due to the coupling 
with 7 . The correction is comparable to the unperturbed effect when ac ~ 5^. Recall that a 
and c are the diagonal transport coefficients and b are the off-diagonal coefficients, relating 
(3 and 7 . Therefore this is a condition on how close the transport coefficients are to a block 
diagonal form, herefore, we have learned that 

The Onsager-Machlup functional, restricted to the quantities [3 which have no external in¬ 
teractions, approximates the probability of a particular relation to a stationary state up to 
corrections which are suppressed by the square of those transport coefficients which mix the 
quantities with and without external interactions. 


V. EXAMPLE: NONLINEAR DIFFUSION 

As an example of the minimum dissipation principle at work, consider diffusion in an 
anisotropic box extending from a; = 0 to x = L in the x direction and in some small 
interval along the y direction. On both boundaries of the x coordinate, couple the system 
to a reservoir. Let the temperature of the reservoir on the right be A units hotter than 
that of the reservoir on the left. As a result there will be a flux of heat from right to left 
across our anisotropic medium. For simplicity, let A be much smaller than either of the two 
temperatures. On the y boundaries place an insulator. This conhguration is depicted in 
FIG.E 

This two-dimensional heat conduction problem involves two thermodynamic forces, the gra- 
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Insulating Walls 



Anisotropic 

Medium 



FIG. 1: A rectangular anisotropic medium is placed between two reservoirs at different but 
similar temperatures. The other two walls are insulating. Heat flows between the two reservoirs, 
establishing a temperature gradient in the x direction. The transport coefficients are nondiagonal 
in the x — y basis, and so there is also a gradient in the y direction. Due to the nonlinearity of 
the flux-force relation in this example, the angle of the isothermal lines depends on the difference 
between the temperatures of the reservoirs. 


dients 


^~dxT ^~dyT 


(27) 


and two currents, which are the heat fluxes Jx and Jy in the x and y directions. Consider 
the asymmetric nonlinear transport equations 


Jx — ^XxT'^Xx Jy — ^yxLT^Xl -|- XyyT'^Xy. 


yyj. ^^y. 


(28) 


Hence, the transport coefficients read 


Qxx ^xxT , Qxy 0 Qyx ^yxLT Xx-, Qyy XyyT 


(29) 


Here the nonlinearity comes from the off-diagonal Xyx term, which creates a temperature 
flux along the y direction as heat flows between the reservoirs. Let Xyx be large, so that the 
nonlinear effects will not be drown out at big T. As the temperature difference between the 
reservoirs is small with respect to the temperature, there will be a stationary state when the 
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temperature gradient in the x direction is constant 




Now the heat flux in the y direction is 


A2\ 

j _ ^ ^yx . r^2Y 
^11 — I-^ “r '^yyJ- -^y' 


LT 


( 30 ) 


(31) 


As we are interested in boundary conditions in which no heat can escape from the boundaries 
in the y directions, one may impose Jy = ^ and so 


A^\y^ _ 1 dT 

" LT^Xyy T2 dy 

and so the temperature gradient along the y direction is 


(32) 


dT 

dy 




yx 


LTX 


(33) 


yy 


As there is no current along the y direction, and as the current in the x direction is constant, 
the divergence of the current vanishes and so this is a steady state. Explicitly one finds 


dT _ „ 5 (AX,,\ 

dt dx\ L ) 



0 


(34) 


While the forces in both the x and y directions are non-vanishing at the steady state, the 
force in the direction 

y = Ty- A^x (35) 

is zero. The temperature in this direction is roughly constant. Notice however, that unlike 
the linear situation, the direction y itself depends on the boundary conditions A. In this 
case we can factorize the problem into two one-dimensional problems, one in the x direction 
and one in the y direction. The y problem is just that of Onsager and Machlup, but to the 
aforementioned effects that come from interactions with the other system if one perturbs 
too far away from the stationary state. Therefore, for sufficiently small perturbations of the 
y profile, the one-dimensional relaxation obeys the minimum dissipation principle. More 
precisely, the relaxation of a perturbation which depends only on y extremizes M varied 
with respect to temperature fluxes along y. In this case we were very fortunate, because 
the two-dimensional problem globally factorized into two one-dimensional problems. This is 
because the integral curves of the y vector form a linear subspace as y is the same direction 
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everywhere. In general, the direction y will depend on the x coordinate. In this case one may 
still factorize problem, but the integral curves whose tangent vectors are y will no longer 
be straight lines as the y direction will be position-dependent. Thus this factorization only 
applies locally, but this is sufficient as we have assumed local equilibrium. Things become 
even more complicated above two dimensions, if the number of directions along which one 
varies is at least equal to two. In this case the vanishing fluxes define a set of vectors at each 
point, but this set of vectors does not necessarily form the tangent space to any foliation 
of our space and so, if some integrability condition is not satisfied, the problem cannot be 
globally factorized into a near-equilibrium problem and a constant problem. However locally 
this factorization is always possible and the principle of minimal dissipation follows. 

VI. CONCLUSION 

We have generalized Onsager’s minimum dissipation principle to relaxations to steady 
states which are only locally in equilibrium. We have used Onsager and Machlup’s stochas¬ 
tic method to demonstrate this generalization of the principle, but we were only able to 
demonstrate it for variations with respect to a subset of the variables for which fluxes are 
not driven by the boundary conditions. However it appears that the full reciprocal relations 
are consistent with the vanishing of variations with respect to all of the fluxes. It should be 
noted that we have not determined quantitatively to which order in the size of the variation 
the principle holds. 

Our generalization differs from that of Ref. [21], in that we allow nonlinear flux-force rela¬ 
tions while they allow microscopic irreversibility. It would be interesting to see if one may 
formulate a generalization which incorporates both theories. This is impeded by the fact 
that our stochastic approach is difficult to generalize to the irreversible case, while their 
Hamilton-Jacobi approach assumes linearity in a number of places, such as the constancy of 
their diffusion matrices and their quadratic ansatz for the Lagrangian density. The next step 
is to incorporate this principle into the geometrical thermodynamic field theory of Ref. [22] . 
To do this, it would be useful to relate M to the length of the system’s trajectory in some 
space, so that the minimum dissipation principle becomes the shortest path. 
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